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We investigate theoretically skyrmion magnonic crystal, i.e., the dynamics of the magnetization in 
a chain of the ferromagnetic nanodots being in skyrmion magnetic configuration. We show that col¬ 
lective excitations are possible to be observed in the structure. We present the dispersion relation of 
the coupled skyrmions. It exhibit a periodical property in dependence on wave vector, characteristic 
feature of the band structure in magnonic crystals. Spatial analysis of the magnetization amplitude 
associated with the magnonic bands confirms type of the excited modes, as breathing and clockwise 
gyrotropic dynamical skyrmions. These high and low frequency excitations, propagate with negative 
and positive group velocity, respectively, and can be explored for study fundamental properties and 
technological applications in spintronics and magnonics. 

PACS numbers: 75.70.Kw, 75.30.Ds, 75.40.Gb, 75.75.-c, 75.78.Fg 


Properties of the magnetic structure are determined 
by set of numerous parameters, like a saturation magne¬ 
tization, exchange constant, anisotropy, damping factor 
or geometry. Because of this large degree of freedom in 
manipulating of the magnetic structures, the physics of 
magnetic materials is extremely vast, rich in phenomena 
and useful for application in microwave technology, infor¬ 
mation storing and processing. Another parameter which 
has important impact on properties of the magnetization 
arrangement and its dynamics is Dzyaloshinskii-Moriya 
interaction (DMI). The DMI is also called an antisym¬ 
metric exchange interaction and it appears as a result of 
spin-orbit coupling in non-centrosymmetric magnets or 
it comes from spin-orbit interactions at the interface of 
ultrathin magnetic films 

Although influence of the DMI on magnetization tex¬ 
ture at room temperature has not yet been demonstrated, 
its influence on dynamical properties is significant. Re¬ 
cently spin waves (SWs) propagating in thin films have 
been investigated analytically^ numerically^ and also 
experimentally^ showing the great importance of DMI 
to nonreciprocal character of SW propagation and elec¬ 
tric field controllable phase shift of SWs. In addition, 
DMI can destroy the degeneracy of chitrality of the vor¬ 
tex stated 

At low fields and temperatures presence of DMI leads 
to appearing of the exotic spin textures. As result of si¬ 
multaneous contribution with exchange interactions, var¬ 
ious arrangements of magnetization can spontaneously 
realize, e.g., helical^ or conical^ spin textures. An inter¬ 
esting phenomena recently discovered in magnetic struc¬ 
tures is a skyrmion^ It was shown with Monte Carlo 
and micromagnetic simulations that single skyrmion can 
nucleate in the isolated disbii 13 or skyrmionic spin tex¬ 
ture can appear in ultrathin film, i.e., skyrmionic crystal 
(SkX). The first experimental observation of SkX was 
made with the use of neutron scattering^ after the theo¬ 
retical predictions^ Analytic, numerical and experimen¬ 
tal investigations are conducted at the present moment 
in order to understand the formation and properties of 
SkX lattices. Material properties are optimized for stabi¬ 


lization of SkX phase at large range of temperature and 
magnetic field. Positive influence of anisotropy on sta¬ 
bilization was also showed^ It was found that dimen¬ 
sion plays important role in stabilization of the SkX: in 
two-dimensional ultrathin films can exhibit the SkX tex¬ 
ture at wide ranges of temperatures and fields showing 
even near room-temperature formation of the skyrmions 
arrayJi Structuralized and confinement structures are 
also studied^ they not only permit on designing man¬ 
ually SkX lattices, but extends the temperature stabil¬ 
ity range. Since the discovery of skyrmions in magnetic 
structures, they are considered for potential application 
as a part of the memory devices^ magnetic racetracks^ 
and logic devices *21 This is due to their high mobility 
under applied currents. The current needed to control 
skyrmion can be few orders of magnitude lower than for 
control of the domain wall*22 



FIG. 1. (Color online) A structure under consideration. The 
bias magnetic field B e xt is out-of-plane plane and directed 
along the z axis, (a) The 3D plot shows the magnetic config¬ 
uration in a single dot (left image) and the m z (to) component 
of the magnetization (right), (b) The m z (to) component of 
the magnetization is plotted as function of x and y for the 
array of 20 nanodisks with 30 nm diameter. 

The dynamical excitations of skyrmions and SkX 
started to be investigated, too. They can be induced 
with external field, spin waves, electric current or tem¬ 
perature. An examples of conducted studies is interac- 
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tion with magnetic droplet^ and SWs scattering on sin¬ 
gle skyrmion. 23 The appearance of the resonant magne¬ 
tization oscillations in a single skyrmion under influence 
of the external dynamic magnetic field was presented, 
i.e., breathing f —i — clockwise and counter-clockwise 
modes i 24 i 25 Also the melting of SkX lattice under ex¬ 
ternal dynamic magnetic field was shown^ The inter¬ 
action of exotic spin texture with SWs are also studied 
since it was suggested that magnons can be used to con¬ 
trol the motion of skyrmions. 27 The dynamics of coupled 
skyrmions^ and collective excitatior*22 were considered 
as well. The spintronic application of skyrmion-based 
spin transfer nano-oscillator (STNO) was proposed^. 

In this Letter, we numerically study collective 
skyrmion dynamics in the chain of the magnetic nanodots 
with skyrmion texture. Investigated system is analogi¬ 
cal to one presented in Ref. [3l|, where collective excita¬ 
tions of the vortex states were studied. However, in our 
study the dots are in skyrmion state and have almost of 
two orders smaller diameter. We show, that collective 
dynamical skyrmions appear in this magnonic crystals 
due to dipole coupling between nanodots. Moreover, we 
demonstrate, that the dynamic signal can be transmit¬ 
ted through the chain by different bands, which are con¬ 
nected with basic excitations of the isolated skyrmion: 
breathing and gyrotropic type of modes. These high and 
low frequency bands are characterized by positive and 
negative group velocities. 
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FIG. 2. (Color online) Dynamical components of the magne¬ 
tization 5m x , Srriy and 5m z of the clockwise, 0.98 GHz (a) 
and breathing, 13.24 GHz (b) modes obtained using FDTD. 
In the left (right) column the real (imaginary) part of 5m is 
shown. 


We have performed finite-difference time-domain 
(FDTD) simulations with mumax 3 solver— and we have 
used uniformly discretized grid with the size of the cell 
0.6x0.6x1.0 nm 3 . In mumax 3 , the time derivative of M 
is defined as the torque r: dMWd _ r> p as i c micro- 
magnetic properties are described with Landau-Lifshitz 


torque: 


r = 7 1 (M x B eff + a (M x (M x B eff ))), (1) 

1 + or 

where a is the dimensionless damping parameter and 
B e ff is the effective magnetic induction field, which 
can have contributions from: external magnetic field, 
magnetostatic demagnetizing field, Heisenberg exchange 
field, Dzyaloshinskii-Moriya exchange field, magneto¬ 
crystalline anisotropy field and thermal field. For pur¬ 
pose of our study, following contributions to the effective 
field B e ff are taken into account: 

Beff = B ex t T Bdemag T B ex ch T Bf)M T B an f s , (2) 

where B ex f, Bdemag? B ex ch? Bdm, B an f s are external 
magnetic field, demagnetizing, exchange, Dzyaloshinskii- 
Moriya and anisotropy fields, respectively^ 
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FIG. 3. (Color online) The |2DFT| function showing (a) 
the dispersion relation of the collective clockwise gyrotropic 
skyrmion excitation and (b) the dispersion relation of the col¬ 
lective breathing skyrmion excitation. The black squares in¬ 
dicate wavenumbers that are investigated in manuscript. The 
|2DFT| is normalized to maximum value for each graph seper- 
ately. The siganl in (a) is ~ 10 times weaker then in (b). 


We have used material parameters of the system 
with perpendicular anisotropy that correspond closely to 
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Pt/Co structure and have been used in the Refs. || and 
El, apart from diameter of the nanodisks and damping 
factor. These two parameters were changed for enhance¬ 
ment of the dispersive character of collective skyrmion 
excitation and better resolution. The parameters that 
were used are: thickness of the dots d = 1 nm, magne¬ 
tization saturation Ms = 10 6 A/m, exchange constant 
A ex = 1.5 x 10 -11 J/m, uniaxial anisotropy constant 
K u i = 10 6 J/m 3 , damping constant a = 0.0001 and DMI 
constant D xn d = 0.003 J/m 2 . We keep in calculations the 
out of plane external magnetic field fixed to B z = 0.1 T. 
The micromagnetic simulations were performed in two 
stages. In the first step the magnetization configuration 
of initial state was postulated with the Neel skyrmion 
in the center of each dot. The system was relaxed to 
achieve energy minimum. 33 Due to that procedure the 
stabilized magnetization vector distribution ( m x (x , y, to), 
m y (x, y, to), m z (x,y,t d )) was obtained. 
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FIG. 4. 

(Color online) ( 

a) Real part of Sm z , 

Re(Jm z ) 

of the 


collective clockwise gyrotropic mode with wavenumber k = 
0.6 X A (b) Ke(Sm z ) of the collective breathing mode with 
wavenumber k — 0.4 x -. 


In the second stage the simulations of the magnetiza¬ 
tion dynamics were performed. The magnetization dis¬ 
tribution obtained in the first stage were influenced by 
an alternating external magnetic field hdyn that was di¬ 
rected either out of plane (z) or in-plane (along x axis). 
This dynamic field h dyn = h amp sinc[27r/ s i gna i(t - 4hift)] 
is characterized with cutoff frequency /signal = 16 GHz 
and is shifted in time by Ghift = 10 -9 s. The ampli¬ 
tude paramter of the sine signal, h amp was choosen so 
that the maximum absolute value of dynamic magnetic 
field was 0.01. The excitation field h dyn is uniform and 
limited to the dot. In dynamic simulations, we used the 
Dormand-Prince method with adaptive time step control 
for advancing the Landau-Lifshitz equation. The maxi¬ 
mum time step was defined to 10“ 12 s and minimum time 
step to 10“ 15 s. The total time of the simulations was 
limited to 3x 10“ 7 s and dynamic magnetization compo¬ 
nents were saved with sampling interval 3 x 10“ 11 s. Such 
specification allows for calculating Fourier transform to 
frequency domain with resolution below 10 MHz within 
the range slightly exceeding (0,/ signa i). 

The first structure under investigation is isolated disk 
with thickness d = 1 nm and diameter 2 R = 30 nm. 
Dynamics of skyrmions in isolated disk was already in¬ 
vestigated and our simulations are in agreement with 


these resultsThe ground state of the isolated disk is 
presented in Fig. [lja). We find two intensive skyrmion 
modes below cutoff frequency of the signal: mew clock¬ 
wise gyrotropic mode at 0.98 GHz (excited with in-plane 
external magnetic field) and me r breathing mode at 
13.24 GHz (excited with out-of-plane external magnetic 
field) ^ The discrete-time Fourier transform (DTFT) 
was performed for magnetization components, separately 
at each point of the 2D FDTD grid and value of DTFT at 
frequency of interest was plotted as function of x and y. 
The function obtained in this way shows dynamic ampli¬ 
tude distribution (its real and imaginary parts, Re(Jra) 
and Im(Jra), respectively) as presented in Fig. O The 
character of the breathing mode is clear to interpret, 
since it shows that largest changes of the magnetization 
z component in form of a ring shape (Fig.[2jb)), meaning 
that skyrmion is dynamically expanding and compress¬ 
ing. The character of the clockwise mode is more difficult 
to intuitively interpret from profiles. However, it can be 
visualized dynamically according to the formul a 35 ! 36 : 

mi(x , y , t ) = rrii(x, y , t 0 ) + cos(2M)Re(Jra*)+ 

sin(2M)Im(Jra^), 

where index i indicates the x, y and z components. The 
animation showing time evolution of all three modes can 
be found in the supplementary material^ 
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FIG. 5. (Color online) The m z (t\) — m z (to) of collective 
breathing mode excited with the frequency /signal = 12.984 
GHz. 

By extending system to the chain of closely spaced 20 
nano disks, we create the ID magnonic crystal based on 
skyrmions. In Fig. QJb) the distribution of m z (to) in the 
array consisting of 20 nanodisk of 30 nm diameter and 1 
nm thickness is presented. The periodic boundary condi¬ 
tions are additionally implemented within mumax 3 with 
finite number of the repetitions (15 in our simulations) 
along the x direction to serve for minimization effects of 
the demagnetizing fields from the external edges of the 
arrays 

In order to obtain dispersion relation of dynamical 
skyrmions a 2D Fourier transform (2DFT) (from space 
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and time domain (x, t) to wavenumber and frequency do¬ 
main (fc, /)) is performed on magnetization z component, 
m z (x,y,t) at value of y corresponding to center of the 
dots. The calculated dispersion relation of dynamical 
skyrmions in the array of nano dots is shown in Fig. [3j 
This dispersion was obtained with out-of-plane excita¬ 
tion field hdyn, applied locally to the disk in center of 
the chain (disk no. 9). Two magnonic bands are found 
in frequencies below /signal? both have clear periodicity 
(with the period G ) in dependence on the wavenumber 
and Brillouin zone (BZ) boundary at k = it/ a ~ 0.1 
nm" 1 . The finite size of the structure used in simula¬ 
tions is seen as discrete states in the dispersion relation. 
This discrete character might be eliminated by increas¬ 
ing the number of states (increasing number of disks), 
” smearing” of discrate states (increasing damping factor 
value) or by influencing the strenght of interactions and 
distance between discrate states (by changing geometry 
of the structure). 

Frequencies of magnonic bands in Fig. [3] (0.49-0.99 
GHz and 11.88-13.68 GHz) correspond to frequencies of 
the modes in the isolated dot. Thus, it suggests that 
these magnonic bands can be related to clockwise gy- 
rotropic mode and breathing mode branch, respectively. 
In order to determine character of the modes and to visu¬ 
alize the spatial distribution of the dynamic magnetiza¬ 
tion component 5m z , we follow the procedure described 
in Ref. (38j| . We perform 2DFT on m z for each value 
of y seperately. In the next steps we define a ID vec¬ 
tor consisting of values of obtained 2DFT at choosen 
frequency, we choose the wavenumber fc, keep the val¬ 
ues of 2DFT at all wavenumbers with the value differing 
by reciprocal lattice vector G = 2ir/a\ k db nG, where 
n is integer and a lattice constant (in our case a = 30 
nm) and set the zero at other wavenumbers. Finally we 
perform inverse Fourier transformation, resulting in dy¬ 
namic magnetization component for mode with specific 
frequency and wavenumber Sm z . The results are pre¬ 
sented in Fig. 2] (a) and (b). The spatial distribution of 
Re(^m^) in magnonic crystal clearly exhibit similarities 
with R e(Sm z ) and lm(Sm z ) of isolated dots (Fig. [2] (a) 
and (b)). 

The band width of the breathing mode branch is more 
than three times wider than the width of the clock¬ 
wise gyrotropic mode branch. This points at the much 
stronger coupling between breathing oscillations. More¬ 
over, the group velocities have opposite signs in these two 
branches. Those features point at the different coupling 
mechanisms responsible for propagation in each case. We 
showed, that the breathing mode has dominating oscilla¬ 
tion of the m z component of the magnetization (see Fig. 
121(b)). Thus, the coupling between skyrmion excitations 
in nanodots is realized mainly via out-of-plane compo¬ 
nents of M. The lower frequency of collective skyrmion is 
expected for antiphase oscillation of m z in the neighbor 
nanodots (assuming prevailing role of m z in coupling), 
analogically to a chain of parallel magnetic dipoles ori¬ 
ented perpendicularly to the chain axis, where the state 


of lower energy is for antiferromagnetic arrangement of 
dipoles. In magnonic crystals, such oscillations are real¬ 
ized at the BZ boundary (where wavelength of the exci¬ 
tation is equal 2a), while in the BZ center all nanodots 
oscillate in-phase. This explain the frequency decrease 
of the breathing mode with increasing wavenumber from 
zero to BZ boundary, i.e., negative group velocity. For 
the gyrotropic mode the interaction is governed by the ef¬ 
fective net magnetization induced by the core shift, sim¬ 
ilar to the gyrotropic modes in nanodots in the vortex 
stated 

Another interesting point is that in the case of iso¬ 
lated disk we have found in simulations only one mode 
(breathing) when the out-of-plane magnetic field was ap¬ 
plied, whereas in case of the array of dots, two branches 
(breathing and clockwise branch) are excited with out- 
of-plane dynamic magnetic field applied to the single dot. 
We suppose that it is due to generation of demagnetizing 
field with three nonzero magnetization components. It 
means that the demagnetizing magnetic field exerted by 
the oscillations in the excited dot couple to the low fre¬ 
quency gyrotropic mode in the near nanodots. Approx¬ 
imate 10 times smaller intensity of the breathing band 
dispersion then the gyrotropic mode band found in cal¬ 
culations (Fig. [3]) confirms this hypothesis. The further 
investigations are necessary to identify if the process of 
coupling is between breathing oscillations and gyrotropic 
modes or only between forced oscillation in single dot and 
gyrotropic motion. 

The interesting property learned from the dispersion 
relation in Fig. [3] is also possibility for dynamic magneti¬ 
zation transmission through the skyrmion array. This is 
due to non-zero group velocity for k =/= mr/a. We excite 
skyrmion mode with sinusoidal microwave magnetic field 
oscillating at /signal = 12.984 GHz (which correspond to 
k = 0.4(^)) and plot m z (ti) — m z (to) at selected time 
t % = 1.5 xlO 9 s in Fig. [5] We apply the external mag¬ 
netic field only at the first disk. The results presented 
in Fig. [5] confirm propagation of the breathing skyrmion 
mode. 

In summary, we have shown that magnonic band struc¬ 
ture based on skyrmion excitations are possible to form 
in array of nanodisks with DMI. The dispersive character 
and collective properties of the excitation were demon¬ 
strated with micromagnetic simulations. Analogically to 
collective vortex excitation, they are potentially applica¬ 
ble for logic device s 39 ! 40 and can be verified experimen¬ 
tally with time-resolved scanning transmission x-ray mi¬ 
croscopy STXM^l. The research opens a new direction 
in research of skyrmions and might be extended to 2D or 
even 3D artificial structures and also any aperiodic struc¬ 
ture in order to obtain desired functionality. Also the 
tailoring of dispersion propertis of dynamical skyrmions 
might be realized by skyrmion phase control in neigh¬ 
bouring dots, e.g., antiparallel aligment. A promissing 
applications of the propesed ID magnonic crystals and 
collective excitations are related to their combining with 
single dot excitation of the skyrmion dynamics by spin 
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